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Abstract. We discuss image formation in gravitational lensing systems using wave 
optics. Applying the Fresnel-Kirchhoff diffraction formula to waves scattered by 
a gravitational potential of a lens object, we demonstrate how images of source 
O ■ objects are obtained directly from wave functions without using a lens equation for 

gravitational lensing. 
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1. Introduction 

Gravitational lensing is one of the prediction of Einstein's general theory of relativity 
and many samples of images caused by gravitational lensing have been obtained 
observationally pp. Light rays obey null geodesies in curved spacetime and they are 
deflected by gravitational potential of lens objects. In weak gravitational field with thin 
lens approximation, a path of a light ray obeys so called lens equation for gravitational 
lensing and many analysis concerning the gravitational lensing effect are carried out 
based on this equation. Especially, we can obtain images of source objects by solving 
the lens equation using a ray tracing method. As a path of light ray is derived as 
the high frequency limit of electromagnetic wave, wave effects of gravitational lensing 
become important when the wavelength is not so much smaller than the size of lens 
objects and in such a situation, we must take into account of wave effects. For example, 
when we consider gravitational wave is scattered by gravitational lens objects, the wave 
effect gives significant impact on the amplification factor of intensity for waves [21 El II]- 
Another example is direct detection of black holes via imaging their shadows [3 [6]. 
The apparent angular size of black hole shadows are so small that their detectability 
depends on angular resolution of telescopes which is determined by diffraction limit of 
image formation system. Thus it is important to investigate wave effects on images for 
successful detection of black hole shadows. 

Although interference and diffraction of waves by gravitational lensing has been 
discussed in connection with amplification of waves, a little was discussed about how 
images by gravitational lensing are obtained based on wave optics. For electromagnetic 
wave, E. Herlt and H. Stephani [7] discussed the position of images by a spherical 
gravitational lens evaluating the Poynting flux of scattered wave at an observer. They 
claimed that there is a disagreement between wave optics and geometrical optics 
concerning the position of double images of a point source. But they have not presented 
complete understanding of image formation. In this paper, we consider image formation 
in gravitational lensing using wave optics and aim to understand how images by 
gravitational lensing are obtained in terms of waves. For this purpose, we adopt the 
diffraction theory of image formation in wave optics [8] , which explains image formation 
in optical systems in terms of diffraction of waves. 

2. Wave optics in gravitational lensing 

We review the basic formalism of gravitational lensing based on wave optics [1] . In this 
paper, we does not consider polarization of waves and treat scalar waves as a model for 
electromagnetic waves. 
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Figure 1. The gravitational lens geometry of the source, the lens and the observer. 
a <C 1 is the deflection angle, tl is the distance from the observer to the lens object 
and rs is the distance from the observer to the source, Tls = rg — tl- 

Let us consider waves propagating under the influence of the gravitational potential 
of a lens object. The background metric is assumed to be 

als 2 = g„ v dx^dx v = -(1 + 2U(r))dt 2 + (1 - 2U(r))dr 2 } (1) 

where U(r) is the gravitational potential of the lens object with the condition \U\ <C 1. 
The scalar wave propagation in this curved spacetime is described by the following wave 
equation: 

d, ( V / =0<r0„$) = 0, (2) 

and for a monochromatic wave with the angular frequency u, 

(V 2 + w 2 )$ = 4w 2 £/(r)$, (3) 

where V 2 is the flat space Laplacian. 

We show the configuration of the gravitational lensing system considering here 
(Figure [Q. The wave is emitted by a point source, scattered by the gravitational 
potential of the lens object and reaches the observer. We assume the wave scattering 
occurs in a small spatial region around the lens object and outside of this region, the 
wave propagates in a flat space. With the assumptions of the eikonal and the thin 
lens approximation, the Fresnel-Kirchhoff diffraction formula provides the following 
amplitude of the wave at the observer |TJ [3] 

$(T7, A) = — f d ^ eX p [icoS(r), £, A)] (4) 
27Tir LS r L J 

where S(ri, £, A) is the effective path length (eikonal) along a path from the source 
position t] to the observer position A via a point £ on the lens plane 

S( V , & A) = [(* - r,) 2 + rls] 1/2 + [(<• - A) 2 + r\] 1/2 - fc) 
(r/ - A) 2 rLrs f£-A r) - A 
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and we assume that |t7 — A| C rg and |£ — A| <C r^. A constant a represents the 
intensity of a point source. The two dimensional gravitational potential is introduced 
by 



V(*)=2/ dzU(tz). (6) 

Then the wave amplitude at the observer can be written as [TJ [2] 

$(77,A) = $ (»7,A)F(77,A) (7) 

where $o is the wave amplitude at the observer in the absence of the gravitational 
potential U : 

$0(77, A) = ^ exp [iuS (r,, A)] , S (rj, A) = (t? ~ A) ' + r 5 . (8) 

5"o(»7, A) is the path length along a straight path from ry to A. The amplification factor 
F is given by the following form of a diffraction integral 

F(r),A) = [ c/ 2 £exp [iuSi(r],€, A)] , (9) 

r L r LS 2m J 

S l(l7 ,«,A) = p(i^-^) 2 -^) (10) 
2r L s \ r L r s J 

where 5*1(77, A) is the Fermat potential along a path from the source position rj to the 
observer position A via a point £ on the lens plane. The first term in S\ is the difference 
of the geometric time delay between a straight path from the source to the observer and 
a deflected path. The second term is the time delay due to the gravitational potential 
of the lens object. Now we introduce the following dimensionless variables: 

x = T , y = — T , d=[l — , w = v-u, v = — 7T ^ 11 

4o r s £ V r s J 4o r LS r L r s t% 

where we choose £0 as 



to = r L 0E, 0E = \^^. (12) 
V r L r s 

M is the mass of the gravitational source, £0 an d @e represent the Einstein radius and 
the Einstein angle, respectively. Using these dimensionless variables, 



F(y, d) = j d 2 xexp iw Q(a; - y - df - ifi(x)j 



$o(2/, d) = — exp 

r S 




«w I - — I y d I + 



it) 



(13) 
(14) 



In the geometrical optics limit w ^> 1, the diffraction integral (flBjl can be evaluated 
around the stationary points of the phase function in the integrand. The stationary 
points are determined by the solution of the following equation: 

x - y - d-V x tp{x) = 0. (15) 

This is the lens equation for gravitational lensing and determines the location of the 
image x for given source position y. As the specific model of gravitational lensing, 
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we consider a point mass as a gravitational source. In this case, the two dimensional 
gravitational potential is 

ip(x) = In \x\ (16) 

and the deflection angle is given by [j] 

4M 

a = |V^| = — . (17) 

For y = d = 0, the solution of the lens equation ffl5l) is 

\x\ = 1 (18) 

and represents the Einstein ring with the apparent angular radius 6e defined by ( TT2j) . We 
show an example of images obtained as solutions of the lens equation (TH)]) (Figure To 
produce these images, we have assumed an extended source with Gaussian distribution 
of intensity. 
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Figure 2. Images of gravitational lensing by a point mass. The source is assumed to 
have the intensity with Gaussian distribution. From the left to the right panels, the 
source positions are y = 0.0, 0.5, 1.0, 1.5. 



The wave property is obtained by evaluating the diffraction integral (Tl3|) . 
point mass lens potential ( fl6l) . the integral can be obtained exactly 

F(y) = e (V2M|y| 2 +lnW2V/4 r A _ l _ w ^ ^ A _ l _ Wj 1; -i W \y\- 

On the observer plane, an interference pattern appears (Figure 3). 
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Figure 3. Amplification factor for w = 10. 
| Using the original variables, the lens equation is ^ A — v ~^' — = 0. 
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Near y = 0, the distance between adjacent fringes of the interference pattern is 



This fringe pattern is interpreted as interference between double images of a point 
source by the gravitational lensing. The question we aim to raise in this paper is how 
the interference pattern on the observer plane is related to the images of gravitational 
lensing in the geometrical optics limit. The wave amplitude on the observer plane does 
not make the image of the source and we have to transform the wave function to extract 
images. To answer this question, we introduce a "telescope" in the gravitational lensing 
system and simulate observation of a star (a point source) using the telescope. With this 
setup, it is possible to understand how images of a source are formed in the framework 
of wave optics. 

3. Image formation in wave optics 

To establish relation between the interference pattern of the wave and the images of the 
source in the gravitational lensing system, we first consider an image formation system 
composed of a single convex lens and review how images of source objects appear in the 
framework of wave optics [8]. 

3.1. Image formation by a convex lens 

Let us $i n (cc) is the incident wave from a point source in front of a thin convex lens 
and $t(^) is the transmitted wave by the lens (Figure HI). They are connected by the 
following relation 



where T(x) is called a lens transformation function. The action of a convex lens is to 
modify the phase of the incident wave. For a point source placed at z = —f (front focal 
point), the incident wave and the transmitted wave are 




(20) 







where we have used r = a// 2 + \x\ 2 ~ / + \x\ 2 /(2f) assuming |je| <C /. Thus, a convex 
lens converts a spherical wave front to a plane wave front. 
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Figure 4. Wave front modification by a convex lens. 

Using this action of a convex lens for the incident wave and the transmitted wave, 
we can demonstrate the image formation by a convex lens in the framework of wave 
optics. Let us consider the configuration of the lens system shown in Figure 




Figure 5. One lens image formation system. 



We assume the distribution of the source field on the object plane z = —a as &q(x ). 
Using the Fresnel-Kirchhoff diffraction formula, the amplitude of the wave in front of 
the lens is given by 



$i(a;) oc / d^x ^ (x )e turi « / d'x <$> (x ) e 



.la:-] — xn I ^ 



where r\ is the path length from a point on the object plane to a point on the lens plane 
and we have assumed \x\ — x \ <«. The amplitude of the wave just behind the lens is 
given by the relation fT2T|) 

$> v {x 1 ) = t L {x 1 )e-^-^f§ l {x 1 ) 

where t^ is the aperture function of the lens defined by ti{x) = 1 for < \x\ < D 
and ti(x) — for D < \x\. D represents a radius of the lens. With the assumption 
\x\ — x\ b, the amplitude of the wave on the z = b plane behind the lens is 



$ 2 (a?) oc / d 2 %i <3>i'(a?i)e 
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For a value of 6 satisfying the following relation (the lens equation for a convex thin 
lens) , 

1 1 _ 1 
the wave amplitude becomes 



(23) 



$ 2 (a;) oc / d 2 Xod 2 Xi& (x )t L (xi) exp 



. :x x 

-IU[ h — 

a 



oc / d 2 x § (x ) 



2J 1 (uD\x/b + x /a\) 



exp 



Xi 
\X 



exp 

12' 



\XqY_ 
2a 



IUO- 



2a 



(24) 



uD\x/b + £Co/a| 

For cuD 00 limit, the Bessel function in (T2"4"j) becomes the delta function and we 
obtains the following wave amplitude on z = b: 



$ 2 (a;) oc / d 2 x ® (x ) x 5 2 



* o[ -b X 



(25) 



~x x~ 
. a b. 

Thus, a magnified image of the source field appears on the z = b plane. This reproduces 
the result of image formation in geometric optics; we have shown that an inverted 
images with magnification b/a of a source object appears on z = b satisfying the lens 
equation (1231 . 

If we do not take uD — > 00 limit, due to the diffraction effect, an image of a point 
source has finite size on the image plane called the Airy disk [Sj. Its size is given by 



(26) 



This value determines the resolving power of image formation system. For two point 
sources at x = — d/2, d/2, their separation on the image plane is bd/a. To resolve 
them, their separation must be larger than the size of the Airy disk: 



d -> A -B 



(27) 



The lefthand side of this inequality is the angular separation of the sources and 9 
determines the resolving power of the image formation system. 



3.2. Image formation in gravitational lens system 

As we have observed that a convex lens can be a device for image formation in 
wave optics, we combine it with a gravitational lensing system and obtain images by 
gravitational lensing. We consider a configuration of the gravitational lens system shown 
in Figure |6]and examine how the images of the source object appear using wave optics. 
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Figure 6. Configuration of a gravitational lens with a convex lens system. Thin 
orange lines represent paths that contribute to diffraction integrals. 



As the source object, we assume a point source of wave. The amplitude of the wave just 
in front of a convex lens is 

$ L (y,d') = %(y,d')F(y,d'). (28) 

This equation is the same as (j7|). After passing through the convex lens, the wave 
amplitude on the image plane z 2 is given by 



$/(77,A) 



d 2 A'^ L (rj, A')exp 



A'|<A 



. V . 



exp 



— (A - A) 2 

2z 2 1 ; 



(29) 



where Ao denotes the aperture of the convex lens. Using dimensionless variables, the 
wave amplitude on the image plane is 



$,(y,d) = - / d 2 d'F(y + d' 

r S J\d'\<do 



(30) 



x exp 



iw 



fr L r s 
\2r 



LS 



1 1 1 

rs z 2 J 



d' 2 -[y + 



r L r s 
r L sZ 2 



d) d' + 



1 ( r L s 



r L r LS z 2 



If we choose the location of the image plane z 2 to satisfy the following "lens equation" 
for a convex lens, 

1 1 1 

rs z 2 /' 

then the wave amplitude on the image plane becomes 



r S J\d'\<d 



d 2 d'F(y + d') exp 



-iw [ y + - — - ol d 



x e <(™/2Mi/,d) (31) 



where g = (r LS /r L )y 2 + {r L r s /r LS z 2 )d 2 . Thus the wave amplitude on the image plane 
is the Fourier transform of the amplification factor F which gives the interference fringe 
pattern. Under the geometrical optics limit w 3> 1, x integral in the amplification factor 
ffTB"]) can be approximated by the WKB form 

F(y + d!) w Ae iw \-^ x *- y - d ' )2 - i ' ( - x * ) \ 
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(32) 



We have assumed that the aperture of the convex lens is sufficiently smaller than the size 
of the gravitational lensing system and \d'\ < d <^ 1 holds. Then, the wave amplitude 
on the image plane is 

$/(j/,d) oc Ae iw ^ {x *- y)2 ^ (x ^e lW9M/2 



d 2 d' exp 



\d'\<do 

x 2nd; 



— iw x* H a 



Ji(it;|ac* + (3d\d ) 



rLsf 

(3 



d' 



~ w\x* + /3d\d ' " r LS f 
For wdo — > oo limit (large lens aperture limit or high frequency limit), we obtains 

®i(y,d) oc 5 2 



x*(y) h d 



Thsf 



(33) 



(34) 



and the image of the point source appears at the following location on the image plane 
determined by the lens equation (|32|) : 

risf 



d 



x x*[y). 



(35) 



(|34|) and (|35|) reproduce the same result of image formation in the geometrical optics 
(ray tracing) in terms of the wave optics. This is what we aim to clarify in this paper. 
If the lens equation fl32l has multiple solutions x*,j = 1, 2, • ■ -, the wave amplitude on 
the image plane becomes 



_ , . 2J 1 (w\x i J ) + (3d\d ) 

®l(y, d) (X > Aj yrr ! 



(36) 



w \xi 



where Aj are constants. 

As an example of image formation in a gravitational lensing system using wave 
optics, we present the wave optical images of a point source by the gravitational 
lensing of a point mass (Figure 7). They are obtained by Fourier transformation of the 
amplification factor F (equation f l3T|) ) and teh lens equation for gravitational lensing 
( 1X5]) has not been used. This procedure corresponds to image formation by a convex 
lens. These images correspond to images obtained by geometric optics (Figure 2). We 
can observe wave effect in these images. 
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Figure 7. Wave optical images of a point source by the gravitational lensing of a point 
mass. Parameters are w = 40, do = 0.5(aperture of a convex lens), y = 0, 0.5, 1, 1.5. 
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In each images, we can observe concentric interference pattern which is caused by 
finite size of the lens aperture and this is not intrinsic feature of the gravitational 
lensing system. We can also observe radial non-concentric patterns. They are caused 
by interference between double images and represent the intrinsic feature of the 
gravitational lensing system. For y = case which corresponds to the Einstein ring 
in the geometrical optics limit, we can observe a bright spot at the center of the ring, 
which is the result of constructive interference and does not appear in geometric optics. 
For sufficiently large values of wdo, the wave amplitude at the observer coincides with 
the result obtained by geometric optics. 

4. Summary 

We investigated image formation in gravitational lensing system based on wave optics. 
Instead of using a ray tracing method, we obtained images directly from wave functions 
at the observer without using a lens equation of gravitational lensing. For this purpose, 
we introduced a "telescope" with a single convex thin lens, which acts as a Fourier 
transformer for the interference pattern formed at an observer. The analysis in this 
paper relates the wave amplitude and images of the gravitational lensing directly. In 
the geometric optics limit of waves, images by lens systems are obtained by a lens 
equation which determines paths of each light rays. As light rays are trajectories of 
massless test particles (photon), expressing image in terms of wave is to express particle 
motion in terms of waves. 

As an application and extension of analysis presented in this paper, we plan to 
investigate gravitational lensing by a black hole and obtain wave optical images of 
black holes. This subject is related to observation of black hole shadows [5j |6]. As the 
apparent angular sizes of black hole shadows are so small, the diffraction effect on images 
are crucial to resolve black hole shadows in observation using radio interferometer. For 
SgrA*, which is the black hole candidate at Galactic center, the apparent angular size 
of its shadow is estimated to be ~ 30/i arc seconds and this value is the largest among 
black hole candidates. For a sub-mm VLBI with a baseline length D, using (127|) . the 
condition to resolve the shadow becomes 

D > 1000 km, 

and this requirement shows the possibility to detect the black hole shadow of SgrA* using 
the present day technology of VLBI telescope. Thus, analysis of black hole shadows 
based on wave optics is an important task to evaluate detectability of shadows and 
determination of black hole parameters via imaging of black holes. 

The topic of wave optical image formation in black hole spacetimes belongs to a 
classical problem of wave scattering in black hole spacetimes [9J. As is well known, waves 
incident to a rotating black hole are amplified by the superradiance [10] due to dragging 
of spacetimes. This effect enables waves to extract the rotation energy of black holes. 
On the other hand, it is known that particles can also extract the rotation energy of 
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black holes via so called Penrose process. By investigating images of scattered waves by 
a rotating black hole, we expect to find out new aspect or interpretation of phenomena 
associated with superradiance in connection with the Penrose process. 
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